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We study homogeneous and isotropic cosmologies in a Weyl space- 
time. We show that the field equations can be reduced to the Einstein 

O" 1 

^ | equations with a two-fluid source and analyze the qualitative, asymp- 

kty totic behavior of the models. Assuming an interaction of the two flu- 

ids we impose conditions so that the solutions of the corresponding 
j> , dynamical system remain in the physically acceptable phase space. 

| We show that in Weyl integrable spacetime, the corresponding scalar 

field acts as a phantom field and therefore, it may give rise to a late 
accelerated expansion of the Universe. 
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1 Introduction 



The usual approaches for an explanation of the late-time acceleration of 

*_ , 

c3 ■ the universe are characterized by a departure from conventional cosmology. 

The proposed models, either assume the existence of dark energy [H [2] , or 
require a modification of general relativity at cosmological distance scales 
[3J H], (cf. [51 [8] for comprehensive reviews and references). Less 
explored is the idea that the geometry of spacetime is not the so far assumed 
Lorentz geometry (see for example [9]). Due to its simplicity Weyl geometry 
is considered as the most natural candidate for extending the Lorentzian 
structure. 

We recall that a Weyl space is a manifold endowed with a metric g and 
a linear symmetric connection V which are interrelated via 

= -Q^ga/3, (1) 

where the 1-form is customarily called Weyl covariant vector field (see 
the Appendix in [TD] for a detailed exposition of the techniques involved in 



Weyl geometry). We denote by D the Levi-Civita connection of the metric 

9a/3- 

A consistent way to incorporate an arbitrary connection into the dynam- 
ics of a gravity theory is the so-called constrained variational principle [11] . 
Applying this method in the context of Weyl geometry to the Lagrangian 
L = R, one obtains the field equations 



where G^ u ) is the symmetric part of the Einstein tensor (see equations (30) 
and (31) in [H]). If we express the tensors Gu^ and M^ v in terms of the 
quantities formed with the Levi-Civita connection D, the field equations 
become 



In the case of integrable Weyl geometry, i.e., when = d^, the source 
term is that of a massless scalar field. Taking the divergence of (T5]) and 
using the Bianchi identities we conclude that 



In this paper we study Friedmann- Robertson- Walker (FRW) cosmologies 
in a Weyl framework. In Section 2 we explore the field equations derived 
from the Lagrangian L = R + L m , where the matter Lagrangian, L m , is 
chosen so that ordinary matter is described by a perfect fluid. It is shown 
that the presence of the Weyl vector field can be interpreted as a fluid and 
we analyze the asymptotic behavior of the models. Assuming an energy 
exchange between the two fluids we extend previous work [TO]. In Section 
3 we consider a modification of the Einstein-Hilbert Lagrangian, cf. (1171) . 
which may provide a mechanism of accelerating expansion. 

2 Interacting fluids 

In the following we assume an initially expanding FRW universe with expan- 
sion scale factor a (t) and Hubble function H = a/a. We adopt the metric 
and curvature conventions of [12]. An overdot denotes differentiation with 
respect to time t, and units have been chosen so that c = 1 = 8ttG. Ordinary 
matter is described by a perfect fluid with energy-momentum tensor, 




V (M g.) + Q^Q U + (y a Q a - Q a Q a ) =: M t 





D^Q^ = 0. 



Tfj, v = (P2 + P2) U^Uu + Pig^u, 



(3) 



supplemented with an equation of state p 2 = (72 — 1) Pi- Since for spatially 
homogeneous and isotropic spacetimes there is no preferred direction, 
must be proportional to the fluid velocity w M , i.e., 

Q"=:gu", Q 2 = Q^ = -q 2 . 

(In vacuum, we have to make the assumption that is hypersurface or- 
thogonal, i.e. it is proportional to the unit timelike vector field which is 
orthogonal to the homogeneous hypersurfaces). Formally the right-hand 
side of (El) can be rewritten as 



I (q^Qu - ^Q 2 g^ = (pi + Pi) u^u u + pxg^, 



(4) 



with 

p\=P\ = -^q 2 , (5) 

i.e., the equation of state of the q— fluid corresponds to stiff matter. There- 
fore we are dealing with a two-fluid model with total energy density and 
pressure given by 

P = pl+P2, P = Pl+P2, (6) 



respectively, where 



Pi = Pi, P2 = (72 - 1)P2, (7) 



i.e., 71 = 2 and 72 < 71. 

The field equations are the Friedmann equation 

H 2 + - 2 =\(pi+P2), (8) 
and the Raychaudhuri equation 

H = -H 2 - \ [(3 7 i - 2) p x + (3 72 - 2) p 2 ) . (9) 

The Bianchi identities imply that the total energy-momentum tensor is 
conserved, so that an interaction between the two fluids is induced. It is 
necessary to make an assumption about the interaction between the two 
fluids (cf [13J), otherwise the field equations constitute an underdetermined 
system of differential equations. The simplest assumption is that the energy- 
momentum of each fluid is separately conserved, so that the two fluids do 
not interact and the densities decay independently, 



Pi = -37i Hp u p 2 = -3<y 2 Hp2- 



(10) 



The dynamical system (j51)- (ITUl) was analyzed in [10J. It was found that 
in expanding models the "real" fluid always dominates at late times and 
therefore the contribution of the Weyl fluid to the total energy-momentum 
tensor is important only at early times. The purpose of this section is to 
weaken the requirement of separate conservation of the two fluids. 

In many cosmological situations the transfer of energy between two fluids 
is important, so one may assume that the two fluids exchange energy. The 
following simple model was proposed by Barrow and Clifton (see [H] for 
motivation and further examples), 

Pi = -37iFpi - f3Hp 1 + aHp 2 , (11) 

p 2 = -3 l2 Hp 2 + 0H Pl - aHp 2 , (12) 

where a and /3 are constants so that the total energy is conserved (see also 
[T5] for a singularity analysis of the master equation derived in |14j). 

Remark 1 In the case of separately conserved fluids, equations / TJOj) imply 
that the sets pi = and p 2 = 0, are invariant sets for the dynamical system 
and by standard arguments, if pi > 0, % — 1, 2, for some initial time to, then 
Pi(t) > throughout the solution. This fact can be made more transparent 
by the following argument. Assuming that 71 > 72, we define the transition 
variable x6 [—1, 1] 

P2 ~ Pi , 1Q x 

P2 + Pi 

which describes which fluid is dominant dynamically /Iff) /. Applying the 
conservation equation to p\ and p 2 , one obtains the evolution equation of 
the variable x, 

X = -#(71-72) (1-X 2 ), (14) 

which implies that the sets \ = ±1 are invariant under the flow of the 
dynamical system. Furthermore, the transition variable x is bounded, that 
is, if initially % G [—1, 1], it remains in that interval for all t. However, the 
choice Ul\) and has the peculiarity that the sets p\ = and p 2 = 0, 

are no longer invariant sets for the dynamical system and therefore the sign 
of the functions p\ and p 2 is not conserved. This is also reflected to the fact 
that the transition variable x n # longer satisfies (Lfy and eventually escapes 
outside the interval [—1,1], thus exhibiting unphysical behavior. 



In order to circumvent these difficulties, one has to impose further con- 
ditions on a and 0. It turns out that the assumption (3 = —a in fTTT]) and 



f[T2"j) is a sufficient condition for the boundness of the function \. With this 
assumption we adopt the Coley and Wainwright formalism [13], for a gen- 
eral model with two fluids. The state of the system consists of the couple 
(x, fi) , where Q = Qi + ^2 is the total density parameter, Q = p/3H 2 . In 
order to allow for closed models in our analysis, we define the compactified 
density parameter co, (see [T6] ) 

n = (15) 

tarn co 

or 

co = arctan [ — - — | , with — n/2 < co < ir/2. 

We see that co is bounded at the instant of maximum expansion (H = 0) 
and also as p — > 0, in ever-expanding models. Finally, defining a new time 
variable r by 

dr 3 (71 - 72) fp 1 



dt 2 V 3 cos co ' 

one obtains the following dynamical system 

duo 1 

— — = — (0 — y) cos 2co cos W 
dr 2 



= (X«~X Jsinw, (16) 

where the constant b is 

6= 3( 7 i+7 2 )-4 ^ 1 
3 (71 - 72) 

In our case, we always have 71 = 2. The parameter Xa which determines 
the range of x is given by 



M 4 « n ^ , 3 (71 - 72) 
Xa = \ 1 - ^7 r, < a < . 

3 (71 - 72) 4 

The phase space of the two-dimensional system ( [TBI) is the closed rect- 
angle 

D=[-7r/2,7r/2] x [-x«,X«] 

in the w — x plane (see Figure 1). Since < Xa < 1, the rectangle D is 
shrinked compared to the phase space in [10] and [TB] . 

The invariant sets of the system are denoted in the following table. 



X 



o- 



Fo Xo 



k = -1 
O < 1 



-tt/2 



fc = +1 



-tt/4 



o- 



H < 



-Xo 



> 1 



= 1 



fc = -1 
n < i 



tt/4 



v ► 

tt/2 u 



fi 2 < fil 



F+ 0+ 
H > 



Figure 1: The invariant sets and equilibrium points of ([16 
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It is easy to verify that the equilibrium points lie at the intersection of these 
sets and are denoted by F x 2 (expanding or contracting flat model) and 
(expanding or contracting open model). The subscripts indicate which fluid 
dominates. Linearization around the equilibrium points is sufficient for the 
characterization of their stability and the result is the phase portrait is 
shown in Figure 2. 

Regions III and IV correspond to expanding models. The F± is a past 
attractor of all models with ft > 0, i.e., the evolution near the big bang is 
approximated by the flat FRW model where the Weyl fluid dominates. Open 
models expand indefinitely and approach at late time a "scaling solution" 
where the Weyl fluid keeps a small fraction of the total energy density. Flat 
models expand indefinitely and the evolution is approximated by the flat 
FRW universe at late time. In both cases the "real" second fluid dominates 
at late times. On the other hand, any initially expanding closed model in 
region III, however close to _F 2 + , eventually recollapses and the evolution is 




Figure 2: The phase portrait of (Tl6|) with 72 = 1. 

approximated by the flat FRW model where the Weyl fluid dominates. 

We therefore conclude that the Weyl fluid has significant contribution 
only near the cosmological singularities. In expanding models the "real" 
fluid always dominates at late times and therefore the contribution of the 
Weyl fluid to the total energy-momentum tensor is important only at early 
times. 

3 Phantom from pure geometry 

The field equations fl2]) constitute the generalization of the Einstein equa- 
tions in a Weyl spacetime in the sense that they come from the Lagrangian 
L = R. There is however an alternative view, namely that the pair (Q, g) 
which defines the Weyl spacetime also enters into the gravitational theory 
and therefore, the field Q must be contained in the Lagrangian indepen- 
dently from g. In the case of integrable Weyl geometry, i.e. when = <9 M 
where is a scalar field, the pair (0, g^ v ) constitute the set of fundamental 
geometrical variables. We stress that the nature of the scalar field is purely 
geometric. A simple Lagrangian involving the set (</>, g^) is given by 



L = R + £V M Q M + L 



(17) 



where £ is a constant and L m corresponds to the Lagrangian yielding the 
energy-momentum tensor of a perfect fluid. Motivations for considering 
theory (IT?]) can be found in [T7J [18] (see also [19] for a multidimensional ap- 
proach and [20] for an extension of ([TTjl to include an exponential potential 
function of (f). By varying the action corresponding to (TT7l) with respect to 
both g M „ and one obtains 

= (<9 M <R0 - i (<9 a 0«9» ^ + T^, (18) 



and 



where □ is the D'Alembertian operator formed with the Levi-Civita connec- 
tion D. As mentioned above, ordinary matter described by T^ u is a perfect 
fluid with energy density p and pressure p. Setting 

we note that for A < the field equations are formally equivalent to general 
relativity with a massless scalar field coupled to a perfect fluid. In Weyl 
spacetime the scalar field has a geometric nature and no restriction exists 
for the sign of the value of A. For A > 0, the Weyl field <fi plays the role of 
a phantom scalar field and therefore, it may provide a mechanism of late 
time acceleration. Further investigation of this issue is the subject of future 
research. 
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